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Abstract 

Tilings of a quadriculated annulus A are counted according to volume 
(in the formal variable q) and flux (in p). We consider algebraic properties 
of the resulting generating function &a(p,q)- For q = — 1, the non-zero 
roots in p must be roots of unity and for q > 0, real negative. 

1 Introduction 

A quadriculated surface is a juxtaposition of squares along sides forming a surface 
(usually with boundary) such that vertices of squares which lie in the interior of 
the surface belong to precisely four squares and vertices in the boundary belong 
to finitely many squares. A simple example of a quadriculated surface is the 
quadriculated plane, M 2 divided into unit squares with vertices in Z 2 . Compact 
quadriculated surfaces have finitely many squares and will usually be bicolored 
and balanced, i.e., the forming squares are black and white in such a way that 
squares with a common side have opposite colors, and the number of black squares 
equals the number of white squares. 

Let A be a balanced quadriculated annulus: A is topologically S 1 x [0, 1] and 
can be embedded in M 2 . We assume A to be a subset of the plane, allowing free 
use of an orientation of A and of its inner and outer boundaries. 

Let Ta be the set of domino tilings of A. For t e T4, we draw inside each 
domino a domino arrow from the center of the black square to the center of the 
white square. Given £o,*i £ Ta, juxtapose the domino arrows from £1 to the 
reversed domino arrows from to'- this produces a finite number of oriented simple 
closed curves. The flux <f>(ti\ t ) of t\ relative to t counts such curves with weights 
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— 1, or 1 given by the winding of each curve in the annulus A. The word flux 
has been chosen for its (co)homological connotations; unfortunately, this use of 
the term is unrelated to that of Lieb and Loss (0). The flux polynomial of A is 

ter A 

this is a Laurent polynomial in p which is independent of base tiling £ 0) up to 
multiplication by an integer power of p. 

Theorem 1.1 All non-zero roots of the flux polynomial $j[(p) are real negative 
numbers. Furthermore, they are either all equal to —1 or all distinct. 

When the roots of 3>j[(p) are distinct, —1 may or may not be a root. 

A tiling t induces a bijection /3 t from black to white squares. Again fixing an 
arbitrary base tiling to G Xa, let the sign a(t; to) of a tiling t relative to t be the 
sign of the permutation Pt°Pt 1 on the set of white squares. Deift and Tomei ([[[[]) 
pointed out that — 1-counting of tilings on disks can only yield very few values. 

Theorem 1.2 (Deift-Tomei) Let D be a quadriculated disk with a base tiling 
to- Then 

5>(*;*o) 

t&T D 

equals 0, 1 or —1. 

The signed flux polynomial of an annulus A is 

site) = to)**** , 

t&T A 

an expression independent of base tiling up to sign and multiplication by a power 
of p. On quadriculated annuli, —1-counting is also severely restricted. 

Theorem 1.3 All non-zero roots of the signed flux polynomial $^(p) are roots 
of unity. 

Actually, we study a more general polynomial, the q-flux polynomial of A, 

<s> A (p,q) = ?" (t;to y (t;<o) , 

teT A 

which enumerates tilings of A according to volume v and flux. Counting according 
to volume (defined in Section 4), i.e., g-counting of domino tilings, has been 
considered previously: in 0, for example, Elkies, Kuperberg, Larsen and Propp 
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g-count tilings by dominoes of the Aztec diamond and MacMahon (||) already 
g-counted lozenge tilings. The flux polynomials are special cases of the g-flux 
polynomial: $\(p) = &a(p, 1) and = ®a{p, — !)• 

The main ingredient in the proofs is the fact that $a(p, g) has two distinct 
algebraic interpretations. The first one is rather familiar: <&a(p, q) is the determi- 
nant of a Kasteleyn-type matrix M^(p, q) (see Section 4). The second interpreta- 
tion requires the construction of another matrix: for each flux /, the coefficient 
of pf in $a(p,C[) (a Laurent polynomial in q) is the trace of a connection ma- 
trix obtained as follows (see Section 6). Tear A along a cut formed by sides of 
squares (as in Figure ^) giving rise to a periodic track segment: a (topological) 
quadriculated rectangle A with a marked pair of opposite sides, the left and right 
attachments (see Figures |9| and |TD|). Now, a tiling of A almost induces a tiling on 
A: the complication is that the dominoes of the tiling, when drawn so as to fill 
A completely, may trespass the attachments. Our definition of tiling of a track 
segment shall allow dominoes to trespass attachments (again as in Figure ||). We 
list, for a given value of the flux /, all possible ways in which dominoes can tres- 
pass the attachment (the indices of an attachment). Say there are N such ways: 
the N x N connection matrix C&j(q) has for (i, j)-entry the (/-counting of tilings 
of A trespassing left and right attachments in the way prescribed by indices i and 
j. Rather unsurprisingly, the connection matrix depends in an irrelevant fashion 
on the choice of cut. 

Both Kasteleyn and connection matrices are well behaved with respect to 
taking the quadriculated n-cover A n of an annulus A. As we shall see in Section 5, 
the Kasteleyn matrix associated to A n has for determinant a polynomial <$>A n (p, q) 
whose roots, for fixed q, are, up to sign, the n-th power of the roots of $^(p, q) 
for the same fixed q. Thus, we may study the roots of $>a(p, q) by considering 
different n-covers of A. For connection matrices, the relationship is even simpler: 
Ca^j = (Ca,/)" (Section 6). Here, A is a periodic track segment obtained from 
A by a cut and A™ is the juxtaposition of n copies of A, i.e., the periodic track 
segment obtained by cutting A n . 

The connection matrix C&nj has two useful features. First, the computation 
of each entry is a counting of tilings on a disk: indeed, the indices associated to 
left and right attachments, obtained by the row and column of an entry, indicate 
how to prune A at the attachments yielding a disk D so that tilings of A with 
prescribed indices are tilings of the disk D, which are then g-counted. We may 



then apply Theorem |L2| to infer that, in the case of (— recounting, connection 
matrices can have very few values for their entries. By finiteness, there must be 
different integers no and n\ for which, for all possible values of the flux /, we must 
have (C A ,/(-l)) rao = (C A|/ (-l)) ni and hence $ A n (p,-l) = $ Ani (p,-l). Since 
the roots of $A n (p, — 1) are, up to sign, n-th powers of the roots of <& A (p, —1), this 
essentially proves that the non-zero roots of <& A {Pi ~ 1) are roots of unity — this 
is Theorem |1.3| . The second feature of the connection matrix is the simple fact 
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that, by definition, it has only non-negative entries when q > 0. This permits the 
study of its eigenvalues (and thus, its trace) by making use of the Perron theorem 
for non-negative matrices. 

In Section 2, we introduce wall- free annuli, essentially irreducible quadricu- 
lated annuli with respect to the set of tilings. For annuli admitting walls, the 
theorems have more stringent statements, with simple proofs. In this section, 
we introduce most of the notation used in the paper and a convenient discrete 
version of the Gauss-Bonnet theorem. 

In Section 3, we review the basic facts about height functions, algebraic de- 
scriptions of tilings of a simply connected quadriculated region ([0). Tilings of 
an annulus A give rise to periodic tilings on A°° , which are then studied as height 
functions in the band A°°. 

In Section 4, we define Kasteleyn weights on A. This is a modification of 
Kasteleyn's construction (|4j], ||). In more detail, cut A to obtain a periodic 
track segment A and a Kasteleyn triple of matrices: one accounting for interac- 
tion between black and white squares whose common side is not in the cut, and 
two others which account for interaction of squares across the cut. The triple 
is convenient when studying the effect of juxtaposition and closing of track seg- 
ments. The g-flux polynomial $a(p, q) is the determinant of a matrix obtained 
in a simple fashion from the Kasteleyn triple. We build up from three simple 
homomorphisms — flux, volume and sign — on the (additive) group of 1-cycles on 
A. The homomorphisms have simple combinatorial interpretations and incor- 
porate well the properties required for obtaining (p, g)-counting of tilings out 
of a determinant. In Section 5, we show how Kasteleyn triples and the g-flux 
polynomial behave with respect to n-covers of the annulus A. 



Section 6 describes the connection matrix C^j- The proof of Theorem [L3 
ends the section. Most of Section 7 is dedicated to the positivity of certain matrix 
entries which is needed to use Perron's theorem. We finally prove the following 



generalization of Theorem O . 



Theorem 1.4 Let A be a wall-free annulus. Let q > be fixed. Then all non-zero 
roots of the q-flux polynomial &a{p,c[) are distinct, negative numbers. 

There is a natural symbolic dynamical context for the results in this paper 
which shall not be explored in detail (see || for terminology and basic results). 
For the band A°°, let 7a°°j be the set of tilings of A 00 with flux /; Ta°°j is a 
Cantor set under a natural sequence space topology. Translation of tilings by 
one period corresponds to the shift in sequence space and is a homeomorphism 
from Ta°°,/ to itself. With this structure, Ta°°j is a topological Markov chain 
with transitive matrix and therefore topologically mixing and with a dense set 
of periodic tilings. Furthermore, translation preserves a certain class of measures 
constructed taking into account the combinatorics of tilings. 
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It seems natural to ask if similar results can be obtained for tilings by lozenges 
of triangulated annuli (a lozenge is a juxtaposition of two equilateral triangles 
along a common side of length 1). We indicate two significant difficulties. First, 
there is no obvious counterpart to Theorem |1.2| . Besides, the technical require- 
ments to prove eventual compatibility, indeed, the very concept of wall, have not, 
to our knowledge, been sufficiently studied. Another related problem is the study 
of other surfaces such as tori and surfaces of higher genus with boundary. This 
must involve polynomials in more variables to account for the value of the flux 
across independent cuts. 



2 Geometry and topology of quadriculated 
annuli 

In the introduction we already presented the concepts of bicolored and balanced 
quadriculated surfaces. We always consider disks and annuli embedded in M 2 but 
not necessarily in the quadriculated plane. Let T$ be the set of domino tilings of a 
quadriculated surface S. A path v Vi . . . v n in S is a sequence of adjacent vertices. 
A path v . . .v n turns left (resp. right) at Vi if and v i+ i are vertices of 

a square to the left (resp. right) of both edges tw^i and ViV i+ i, otherwise, if 
Vi is an interior point in S, the path follows straight at t>j. More generally, the 
curvature of the path . . . Vi-\ViVi + \ ... at Vi is 1, or —1 if i>j is in the interior 
of S and the path turns left, follows straight or turns right at i>j, respectively; if 
Vi is a boundary vertex, let n be the number of squares in the angle Vi-\ViVi + \: 
the curvature at is ±(2 — n), with the sign depending on the orientation of the 
angle. Clearly, the definition of curvature at interior vertices can be considered a 
special case of the more complicated definition for boundary vertices. A path is 
simple if its vertices are distinct and closed if Vq = v n , and in this case we write 

"mmodn 1 

A cut is a simple path in a quadriculated annulus joining outer to inner bound- 
ary and having exactly one vertex in each boundary component. We consider an 
annulus as a circular juxtaposition of quadriculated disks with specified attach- 
ments. We draw cuts and call the disks between consecutive cuts track segments 
since we imagine the annulus to be a circular track for a toy train. More precisely, 
a track segment A is a quadriculated disk whose boundary is divided counter- 
clockwise into four arcs ao, Po and a>\ (for inner boundary, attachment 0, outer 
boundary and attachment 1). Attachments will be oriented from j3 Q to this 
coincides with the boundary orientation for a% but not for olq. A track segment 
A obtained from a balanced annulus A by a single cut is periodic. 

The shape of a path v . . . v n is the sequence of the n — 1 curvatures of the 
path at each Vi, i = 1, . . . ,n — 1. The juxtaposition of track segments A and A' 
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is natural: attachments at and a f must have the same shape and their edge-by- 
edge identification yields the larger track segment AA' with attachments «o and 
a[. The annulus A is the closing-up clA of the periodic track segment A and 
is obtained by attaching a\ to ao; the resulting cut is said to be induced by A. 
Figure |2] shows the closing-up of a periodic track segment and the closing-up of 
the juxtaposition of two track segments. 


































-> 


«0 







Figure 1: Closing-up of track segments 

The juxtaposition A n of n copies of a periodic track segment A under closing- 
up yields the n-cover A n = cl(A n ) of A. Similarly, the universal cover A°° 
is also denoted by A z , an infinite quadriculated band. The canonical projection 
7r : A 00 — > A and the translation by one period r : A°° — > A°° are defined as usual. 
Notice that r takes neighboring vertices to neighboring vertices, satisfies nor = n 
and, for any p G A°° and any curve 7 from p to r(p), tt o 7 goes around A once 
counterclockwise. The projection 7r induces a labeling of boundary components 
of A°° as inner and outer as well as orientations on both boundary components; 
we continue to call such orientations counterclockwise or clockwise. The notation 
A z suggests that we see this band as a two-sided infinite juxtaposition of equal 
track segments 

A oo = A z = ---A_ 2 A_ 1 A A 1 A 2 -- - . 

The induced cut between Aj and A i+1 will be labelled Notice that r takes 
Aj to A i+ i and to £ i+ 3. 

Let A be a quadriculated annulus, £ be a cut and t G T4 a tiling of Draw 
domino arrows inside each domino in t from black to white squares. The flux 
4>(t,^) of t relative to £ is the number of such arrows crossing £ counterclockwise 
minus the number of arrows crossing £ clockwise. 
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Lemma 2.1 Let A be a periodic track segment, A = cl A, t G and let A°° 

have induced cuts £ n+ i. T7ie value of <p(t; £ n+ i) is the same for all n e Z. 

Proof: The numbers of black and white squares in Ao are equal. Also, the 
numbers of black and white squares in the union X of all dominoes in t with at 
least one square in A are also equal. But x = <p(t; £± ) — is the number 

of white squares in X — Ao minus the number of black squares in X — Ao and 
therefore x = 0. The figure shows how each domino contributes to x. ■ 





■ 













Figure 2: Comparing fluxes across cuts 

There is a natural concept of translation of a tiling tA°° of A°°: r(tA°°) is 
formed by the images under r of the tiles of tA°°- A tiling tA°° is periodic if 
r(t j4 oo) = t A oo. A tiling t& of A induces a tiling 7r -1 (^) of clearly, a tiling 
of is induced by a tiling of A if and only if it is periodic. 

A path is a zig-zag if it turns alternately left and right (Figure |3] shows four zig- 
zags). Certain quadriculated annuli A admit a simple closed zig-zag (. Keeping 
in mind that A is embedded in the plane, ( splits the plane in an inner and an 
outer region and splits in turn A as a disjoint union of connected components. 
Usually these components consist of two annuli, as will be further discussed in 



Proposition |2.5| If, however, the zig-zag hits the boundary of A, the number of 
components may increase but they will always be disks and annuli; the zig-zag 
may also coincide with one of the two boundary components of A. A simple 
closed zig-zag ( in a balanced annulus A is a wall if the numbers of black and 
white squares of A in the bounded component of R 2 — ( are equal. 

Proposition 2.2 Let A be a balanced quadriculated annulus with a wall (. No 
domino in any tiling of A may trespass ( . 



Proof: Let X be the set of squares of A in the bounded component of M 2 — £ 
and Y D X be the union of all dominoes with at least one square in X. Both X 
and Y are balanced but all squares in Y — X are of the same color. Thus Y — X 
is empty. ■ 
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An annulus A with no walls is wall-free. From the proposition above, if a wall 
C splits A as a union of two annuli A% and A 2 then the flux polynomials satisfy 
®a = ®a®a 2 and $ A = $ Ai <S>a 2 . 

If a wall £ does not touch the inner (resp. outer) boundary of A, then ( may be 
shifted inwards (resp. outwards) to obtain another parallel wall C (resp. £°); we 
may of course repeat the procedure as often as necessary to obtain parallel walls 
( ! and (° which touch the inner and outer boundary, respectively (see Figure 
U). This decomposes A as a union of ladders (annuli between consecutive parallel 
walls) and disks such that any tiling of A restricts to a tiling of each ladder and 
disk. 




Figure 3: Walls are zig-zags 

A ladder in a tiling is a sequence of parallel dominoes side by side such that 
two neighboring dominoes always touch along an edge of the longer side, each 
domino in the ladder has two neighbors in it and these two neighbors touch the 
domino at different squares. Ladders are then either periodic (in an annulus) or 
bi-infmite (in a band). Periodic ladders in a tiling are always contained between 
two walls. The tilings in Figure f| (which may be interpreted as either a band or 
a cut annulus) are the four tilings of the surface consisting of two ladders. Notice 
that inner and outer boundaries of the annulus in Figure |] are walls. 

Theorem 2.3 Let A be a quadriculated annulus admitting n + 1 parallel walls, 
n > 0. Then 

$» = C + (p + ^(p) = C-(p- 1)", 

for constants C + and C~ . In particular, all roots of <&\{p) (resp. Q^ip)) are 
equal to —1 (resp. 1). 

Proof: Since tilings of A do not trespass walls, p-counting tilings of A amounts to 
p-counting tilings of connected components of the complement of walls in A (i.e., 
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Figure 4: Ladders 



of disks and ladders) and multiplying the results. If L is a ladder, = p+1 and 
$^ = p— 1, up to sign and multiplication by an integer power of p. The restriction 
of a tiling of A to a connected component which is a disk bears no relationship 
with the flux of the tiling: disks therefore only contribute with multiplicative 
constants. ■ 



This proves Theorems and L~3 and gives a fairly complete description of 



the set of tilings for annuli which admit walls. From now on we may consider 
wall-free annuli but this hypothesis will only be crucial in the final section. 

The total curvature of a closed path is the sum of the curvatures at its vertices. 
Interior vertices on a zig-zag alternate between curvatures 1 and —1 and closed 
zig-zags have total curvature 0. Paths with curvature at all interior vertices are 
a natural quadriculated version of geodesies and many of the results for zig-zags 
hold true for geodesies. In this paper we have no use for this kind of path and 
we drop the subject altogether. 

Let v G dS be a boundary vertex of an oriented quadriculated surface S: 
the curvature of the boundary at v is 2 — n v if there are n v squares adjacent 
to v . This coincides with the previous definition of curvature of the boundary 
component once it is given the orientation which leaves the surface to the left. 
The total boundary curvature of a compact quadriculated surface is the sum of the 
curvatures of each boundary vertex. The next proposition is the quadriculated 
version of the Gauss-Bonnet theorem; notice that the usual 2% is replaced by 4. 

Proposition 2.4 The total curvature of the boundary of an oriented compact 
quadriculated surface S is Ax{S). 

The Euler characteristic of S is x{S) = F — E + V, the number of faces 
(squares) minus the number of edges (sides of squares) plus the number of (interior 
and boundary) vertices. 
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Proof Let V{ be the number of interior vertices. Counting vertices face by face 
and again edge by edge yields 



n 



2E = W l + ^(ra„ + l). 




veds 



These three identities imply 



A(F-E + V) = Y t ( 2 -rh), 



veds 



the desired result. 



Thus, since annuli have Euler characteristic 0, the curvature of the outer 
and inner boundary components add to 0. The curvature of the outer boundary 
component may assume any integer value. 

Proposition 2.5 Let A be a quadriculated annulus embedded in M. 2 . 

(a) If a closed zig-zag ( exists in A then all zig-zags in A are simple, each bound- 
ary component of A has curvature and another maximal zig-zag intersecting 
( joins both boundary components of A. 

(b) If a self-intersecting maximal zig-zag ( exists then the curvatures of the bound- 
ary components are 1 and —1 and the endpoints of( are both on the boundary 
component with positive curvature. 

(c) For arbitrary A, there exists a simple maximal zig-zag joining both boundary 
components. 

Proof: Closed zig-zags have total curvature zero. Also, zig-zags intersect by 
sharing an edge. For a self-intersecting zig-zag £, construct a simple closed path 
( of minimal length by starting at an edge used twice and following the zig- 
zag until returning to that edge. The path ( has total curvature ±1 and from 
Gauss-Bonnet cannot bound a subdisk of A. For the same reason, simple closed 
zig-zags cannot bound subdisks of A. Thus, closed zig-zags are simple and wind 
once around A. 

Let Pi and f3 be the inner and outer boundary components of A. If there 
exists a closed zig-zag £ then, applying the technique of proof of Gauss-Bonnet 
for the region bounded by ( and j3 Q , we see that the curvatures of ( and f3 Q add 
to zero. Since the curvature of ( is zero, the curvature of both j3 Q and Pi are then 
also zero. Similarly, if a self-intersecting zig-zag exists then the curvatures of (3 
and Pi are ±1. 

To finish item (a), take a closed zig-zag (. From an edge of ( draw the other 
possible zig-zag C" 1 : from Gauss-Bonnet, C" 1 cannot intersect ( a second time and 
it must eventually hit both boundary components. 
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For item (b), use Gauss-Bonnet now to show that the self-intersecting zig-zag 
( cannot self-intersect twice and its endpoints must both be on the boundary 
component of positive curvature. 

For item (c), we make use of wiggling paths: a path i> . . .v n is wiggling if its 
curvature is ±1 at every vertex Vi, < i < n. In other words, a path is wiggling 
if it always has squares of the same color to its left (or right) . An arbitrary path 
v o . . . v n may be transformed in a wiggling path having, say, white to its left and 
with the same endpoints by substituting each edge ViV i+ i of the path with black 
to its left by ViVi'Vi»v i+ i, where these are the four vertices of a square. An edge 
ViV i+ i, < i < n — lina wiggling path is zig-zig if the curvatures at and v i+ i 
are equal. Clearly, zig-zags are wiggling paths with no zig-zig edges. 

Consider now the class of wiggling paths of minimal length joining both 
boundary components. Such paths are simple, as Figure should convince the 
reader. If there is a zig-zag in this class, we found the required zig-zag cut; this 
is not always the case so in some cases we have to work more. 




Figure 5: Wiggling paths may be taken simple 

What is true instead is that if the class contains no zig-zag it contains a path 
with a single zig-zig edge e. Indeed, choose in the class a path such that the length 
from the inner boundary to the first zig-zig is maximal. If the path contained 
a second zig-zig of identical curvature, we could construct a shorter wiggling 
path, a contradiction (see Figure ^). If, on the other hand, the path contained 
a second zig-zig of opposite curvature, we could exhibit another wiggling path 
which is either shorter or has the same length but a longer initial segment, again 
a contradiction (see Figure P). Both constructions indicated in the figure are well 
defined in the sense that the existence of the required extra vertices and edges 
follows from the fact that the stretch of path being modified is initially interior. 

We thus have two simple zig-zag segments Q and ( Q , starting respectively from 
fa and (3 Q with a common edge e. Continue ( Q as a zig-zag beyond e until it hits 
the boundary or self-intersects. If it hits the inner boundary we are done. If it hits 
the outer boundary, then ( Q together with a piece of f3 Q bound a quadriculated disk 
into which the zig-zag continuation of Q enters at the edge e. Thus, from Gauss- 
Bonnet, the zig-zag continuation of Q eventually hits (3 without self-intersecting 
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Figure 6: Improving wiggling paths 



or intersecting the extended ( Q : Q is then the required zig-zag. Finally, if ( Q 
self- intersects, we are in the situation of item (b) and f3 Q has curvature +1. The 
zig-zag continuation of Q cannot self-intersect (from (b)) and cannot intersect £ G 
a second time (from Gauss-Bonnet) and thus must eventually hit (3 and we are 
done. ■ 



Corollary 2.6 Let A be a balanced quadriculated annulus and n a positive in- 
teger. Then A admits closed zig-zags if and only if A n does. Furthermore, A is 
wall-free if and only if A n is wall-free. 



Proof: A closed zig-zag ( in A clearly lifts to another in A n . Conversely, the 
projection of a closed zig-zag ( n in A n to A can not self intersect (by item (a) of 
the previous proposition) and is therefore a closed zig-zag ( in A. Also, let d be 
the number of white squares minus the number of black squares in the bounded 
connected component of M. 2 — £: the corresponding difference for ( n is nd. Thus, 
( is a wall if and only if ( n is. ■ 



3 Height functions 



We recall the basic properties of height functions for a connected and simply 
connected quadriculated surface S with boundary (|F0|; see a l so f° r ^ ne mmor 



adjustments required for surfaces with no boundary). Let Vs be the set of vertices 
of squares in S and Vss ^ be the subset of vertices on the boundary of S. Call 
an arbitrary but fixed vertex Vb G Vos the base vertex. For a tiling t G T$, 
construct the height function 9 : Vs — * Z as follows: 



(a) 9{v b ) = 0; 

(b) if there is a white (resp. black) square to the left of the oriented edge v V\ 
not covered by a domino of t then 9{v{) — 6(vq) = 1 (resp. —1). 
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Clearly, these rules are consistent when constructing 9 along the path surround- 
ing a domino. This local consistency, together with the fact that S is connected 
and simply connected, yield uniqueness and global consistency of 9. Notice that 
all height functions coincide on the boundary and that changing the base vertex 
alters 9 by an additive integer constant. Height functions admit the intrinsic 
characterization below, given in strictly local terms (for a proof and some appli- 
cations, see H). 

Proposition 3.1 Let S be a connected and simply connected quadriculated sur- 
face with a base vertex vj, G Vqs- Let 9 be an integer valued function defined on 
Vs- Then 9 is the height function of some tiling t (relative to the base point Vf,) 
if and only if the following conditions hold: 

(a) 9(v b ) = 0; 

(b) if the oriented edge vqVi is in the boundary of S and there is a white (resp. 
black) square to its left then 9{v\) — 9(v ) = 1 (resp. —I); 

(c) if the oriented edge vqV\ is in the interior of S and there is a white (resp. 
black) square to its left then 9{y\) — 6(vq) = 1 or —3 (resp. — 1 or 3). 

Recall that all height functions 9 with base point V\, equal at Vb- From the 
local characterization of height functions, given a point v G Vs at a distance i 
(measured along edges) from Vb, \9(v)\ < 3£. Thus, even for infinite connected and 
simply connected quadriculated surfaces, height functions are locally bounded 
in the sense that given any finite subset X of Vs the restrictions of all height 
functions to X are bounded by a constant Cx- In particular, the maximum or 
minimum of any (finite or infinite) set of height functions is well defined and, 
again from the local characterization, is a height function. The set of height 
functions, and therefore the set of tilings, thus form a lattice. 

Since the quadriculated annulus A is not simply connected, the construction 
of a height function as above is not necessarily globally consistent. In [[|, this 
difficulty is addressed by the construction of height sections. Here, we follow the 
simpler alternative of considering tilings in the band A°°, which correspond to 
height functions from Va°° to Z. For convenience, we assume the base point v b 
to belong to the outer boundary of A°°. 

Lemma 3.2 For a fixed reference vertex v r on the inner boundary of A°° , 9(v r ) = 
40(t; £i) + c for all tilings t of A°° and some integer constant c. 

Proof: If the base vertex Vb and the reference vertex v r are the two ends of £ i 

2 

then 9(v r ) can be computed from 9(vb) = using the definition of height function 
by following £i. More precisely, if Vi and v$ + i are consecutive vertices on £i then 
9(vi+i) — 9{vi) equals 
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• 3 (resp. —3) if a domino crosses £1 counterclockwise (resp. clockwise) at 
the edge ViV i+ i, 

• 1 (resp. —1) if no domino crosses £i at fjfj+i and the square to the left of 
ViVi+i is white (resp. black). 

Let n 3 , n_ 3 , n\ and n_i be the number of edges in £i in each situation above. 
Clearly, 

9(v r ) = -3n_ 3 - n_i + ni + 3n 3 = (ni + n_ 3 - n_i - n 3 ) + 4(n 3 - n_ 3 ). 

By definition, 0(t; £i) = n 3 — n_ 3 . But ni + n_ 3 (resp. n._i + n 3 ) is the number 
of edges in the cut with a white (resp. black) square to the left and therefore 
c = rii + n_ 3 — n_i — n 3 does not depend on the tiling t. 

In the general case when vj, or t> r are not ends of £i, one has to account for 
the variations of 9 along stretches of boundary from ends to or v r which again 
are independent of tiling. ■ 

Notice that a height function 9 corresponding to a periodic tiling of A°° does 
not usually satisfy (6 o t)(v) = 9(v) but instead (9 o r)(v) = 9(v) + d for some 
integer constant c'; c' depends on A only and not on the tiling t. With a slight 
abuse, we still call such height functions periodic. Call 0/ )max (resp. 0/,mm) the 
maximum (resp. minimum) of all height functions with a given flux /. These 
extremal height functions are periodic. Indeed, from the lattice structure on the 
set of height functions and maximality of 9f !iaax , it is easy to see that 

#/,max = max(6 l / imax , 0/, max or- c')% max o r - c, 
where d = 6 l /jmax (r(i; b )). 

The following result is an application of height functions to tilings of annuli. 

Theorem 3.3 Let A be a quadriculated annulus and let / max and / m ; n be the 
maximum and minimum values of the flux among all tilings of A relative to an 
arbitrary cut £. Then there exist cuts £ max and £ m i n such that no tiling of flux 
/max trespasses £ max and no tiling of flux f min trespasses £ m i n - 

Like walls (in Theorem p.3|) , the non-trespassable cuts above allow for a mul- 
tiplicative principle to be applied to the counting of tilings of extremal flux. In- 
deed, let be the number of tilings of the covering space A n with flux k + f m i n , 

k = 0, . . . , / max - / min ; we have afl x _ /miii = (a/!iax-/ m J n and ^ = ^o 1 )™- More 
general relations among the will be presented in Proposition |5.2| . 
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Proof: We only prove the existence of £ max . We are searching for a cut on 
which all height functions with flux / max coincide: this happens if and only if 

#/ ma x,min(w) = #/max,max(«) for v on the cut. Let 

X = {V G V J 4-|% max ,min(w) = #/ max ,max(v)} 

be the set of vertices where all height functions of flux / max coincide. This set 



contains the outer boundary since Vf, belongs to it and, from Lemma |3.2| , also 
contains the inner boundary. Let Y be the connected component of X containing 
the outer boundary of A°° (i.e., the maximal subset of X containing the outer 
boundary and such that if two vertices V\,v 2 G X are adjacent and V\ G Y then 
v 2 G Y): we must prove that Y also contains the inner boundary. Set 



Q( v ) = jXnax.min^), V G Y, 

" l^ max ,minM+4, veA°°-Y. 



The function 9 is a height function: we check the local conditions in Proposition 



3J]. The function 9 clearly satisfies items (a) and (b). As to item (c), if v\ 
and v 2 are adjacent and either both belong to Y or neither belongs to Y then 
9{v x ) - 9(v 2 ) = 0/ max ,min(fi) - 0f w ,mm( v 2)- If v i G Y and v 2 G" Y then v 2 £ X 
and we have, from the adjacency of vi and v 2 , 9 fin ^ max (v 2 ) - 9 fmaxMn (v 2 ) = 4 
and therefore 0(ux) - 9(v 2 ) = 9^^^) - 9 fiam>iaax (v 2 ) . Hence in all cases 
9(vi) —9{v 2 ) coincides with the difference in values at v\ and v 2 of either 0/ max , m in 
or ^/ maxi max, genuine height functions, thus showing that 9 satisfies (c). If Y 
contains the inner boundary then (f)(9,£,) = / max but otherwise 4>(9, £) = / max + 1, 
contradicting the maximality of / max . I 

Notice that / max — / m i n is the number of non-zero roots of Also, from 

the definition of flow relative to a cut, / max — / min is less than or equal to the 
length of the shortest cut in A. 



4 Matrices and flux polynomials 

Kasteleyn (0]) showed how to count tilings of quadriculated disks by computing 
a determinant; this idea has been extended by many authors, in particular by 
producing a matrix whose determinant g-counts tilings of quadriculated disks 
(0). A variation of these ideas, to be detailed in Proposition [0| obtains the 
q-flux polynomial &a(p, <?) for a quadriculated annulus A. The g-flux polynomial 
will extend the flux polynomials defined in the introduction: $^(p) = <&a(p, 1) 
and $^(p) = $ A (p,-l). 

Let A be a periodic track segment (as defined in Section 2), A = cl A and £ the 
induced cut in A. The adjacency graph G a of the quadriculated annulus A is the 



15 



oriented graph whose vertices are squares of A and whose edges connect squares 
sharing a common side, always oriented from black to white. The domino arrows 
of a tiling are edges in Ga- An edge crosses £ counterclockwise (resp. clockwise) 
if it leaves through the attachment ao (resp. a\) and enters through o>i (resp. 
cto). Let C\(A) be the set of formal sums (with integer coefficients) of edges of 
Ga- Also, let Zi(A) C Ci(A) be the additive subgroup of formal sums such that 
at each vertex of Ga the sum of the coefficients of the adjacent edges is zero. For 
example, the difference c = ti — t 2 of two tilings of a quadriculated surface A 
is an element of Zi(A). Following the homological terminology, the elements of 
Zi(A) are 1- cycles. 

There is a natural basis for Zi(A): the hole basis. Recall that A and therefore 
Ga are embedded in M. 2 . The bounded connected components of the complement 
of the graph Ga are the holes of the graph Ga'- each hole gives rise to an element 
of the hole basis, as follows. Many holes of the graph are bounded by four edges, 
the small holes. The basis 1-cycle s» corresponding to a small hole consists of 
those four edges oriented counterclockwise (here a white-to-black edge is identified 
with minus the same edge with the original black-to- white orientation). The 
only remaining hole, the large hole, contains the bounded connected component 
of R 2 — A. To obtain the 1-cycle £ associated to the large hole, add the edges 
contained in the boundary of the closure of the large hole, orienting them so as 
to have the hole to the left. Figure below represents the same basis element £ 
in two ways, and is given to justify the apparently convoluted definition. 
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Figure 7: The large hole 

As usual, a circuit in a graph is a 2-regular subgraph; thus, a circuit is a 
simple closed curve contained in Ga- Oriented circuits are 1-cycles and generate 
Zi(A): each Sj is an oriented circuit and £ is a sum of disjoint oriented circuits. 
The decomposition in the hole basis of a circuit c oriented counterclockwise is 
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very geometric. From the region of the plane surrounded by c, remove Ga to 
obtain a disjoint union of holes hi, . . . , hp: we clearly have c = ^2 hi- 

We define three group homomorphisms by prescribing their values on the hole 
basis. 

Definition 4.1 On the hole basis, the flux <p : Z\{A) — > Z, volume v : Zi(A) — > Z 
and szon a : Z\(A) — > {±1} /iai>e i/ie values 

• <f>( 8i ) = 0, 0(£) = 1; 

• = 1, u(£) = 0; 

• cr(sj) = —1, o~(£) = (—l) k+1 , where the boundary of the large hole has length 
2k. 

Given the example in Figure [7], the concept of length of I may appear con- 
fusing. Consider the boundary of the large hole: this consists of fc S im P ie edges 
which bound the large hole on one side only and of ^double edges which bound 
the large hole on both sides (in Figure 0, /c simp i e = 20 and fcdoubie = 3). Let 

2 ^-simple 2fcdouble- 

The definitions of flux and sign of a tiling given in the introduction are special 
cases of the new definitions, as will be seen in Corollary [4.3| . For the moment, we 
provide geometric meaning to these homomorphisms. 

Proposition 4.2 For any oriented circuit c, 

(a) the flux 0(c) is the element of Hi(A) = Z associated with c, i.e., 0(c) = —1, 
1 or depending whether c surrounds the annulus's hole clockwise, counter- 
clockwise or not at all; 

(b) the volume z/(c) is the number of small holes surrounded by c, with a negative 
sign if c is oriented clockwise; 

(c) the sign a(c) equals (— i) 6 + w + fc + 1 ; where 2k is the length of c and b and w are 
the number of black and white vertices in the interior of c. 

Proof: Let c be an arbitrary oriented circuit. For flux and volume, the propo- 
sition follows directly from the decomposition of c = hi + • • • + hp in the hole 
basis; we give the details for sign. Let b, w and k be as above. By definition, 
cr(c) = Y\ °~(hi) an d 

o-(c) = (_lf+^ % 

where 2fcj is the length of the boundary of h^ Notice that the total number of 
edges on or inside c is E = k + fcj. Also, the number of vertices of Q on or 
inside c is V = 2k + b + w. Since the Euler characteristic of the disk surrounded 
by c is 1, F - E + V = 1 and F = E -V + 1 = 1 - k - b - w + Y^ih, whence 

o-(c) = (_l)Hw+fc+l_ 
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So far, 0, a and v have been defined in Zi(A). However, in the introduction, 
we defined 4>(t;to) and o~(t;to) for a tiling t relative to a base tiling t - We now 
relate this old definition of (f> and a with the new one. Notice that a tiling t does 
not belong to Z±(A) but the difference t — to does. 

Corollary 4.3 Let t be an arbitrary tiling of a quadriculated annulus A with base 
tiling t . Then <f>(t; t ) = <p(t — t ) and a{t; t ) = a(t — t ). 

Proof: This follows from the proposition keeping in mind that the difference 
between two tilings is a sum of disjoint oriented circuits q surrounding tiled 
regions, so that hi = Wi and <r(cj) = (— l) ki+1 . ■ 

In parallel with this corollary, we also set v{t\ to) = u{t — to), the volume of 
a tiling t relative to a base tiling t - Also, the reader may check that 4>(t\to) = 
4>(t,£) — <f>(to,£) for an arbitrary cut £. 

Let A a commutative algebra with unit over C; in this paper, A is typically 
a ring of Laurent polynomials in one or two variables: A q = C[q, q^ 1 ] and A PA = 
C[p, p -1 , q, q' 1 }- As usual, an attribution of invertible elements of A, i.e., elements 
of the multiplicative group A*, to the edges of Ga is a A* -weight in A and turns 
Ga into a weighted graph. A weight u is thus a function from the set of edges to 
A* and may naturally be identified with a homomorphism uo : C\{A) — > A*. In 
particular, for a tiling t, cu(t) is the product of the weights of the edges in t. An 
oriented circuit c is the sum of its edges which are oriented black-to-white minus 
the sum of its remaining edges: u(c) is the product of the weights of the edges in 
the first class divided by the product of the weights of the remaining edges. For 
a periodic track segment A, we are interested in certain weights in Ga for the 
annulus A = cl A. 

Definition 4.4 Let A be a periodic track segment, A = cl A and £ the induced 
cut. A A* q -weight u in A is a Kasteleyn weight if: 

(a) the weight of an edge interior to A belongs to A*, 

(b) the weight of an edge crossing £ counterclockwise belongs to pA*, 

(c) the weight of an edge crossing £ clockwise belongs to p~ l A* q , 

(d) for any oriented circuit c, cu(c) = a{c)p^q u ^ . 

Proposition 4.5 Any periodic track segment A admits Kasteleyn weights. Fur- 
thermore, given a base tiling to of A = c\A, for any tiling t of A we have 
cu(t) = u(to)a(t - £ o )p 0( *~ to V ( *~* o) . 
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Proof: Pick a maximal tree inside Ga and assign arbitrary weights in A* to 
the edges of this tree; in particular, assigning 1 to all edges of the tree will work. 

Given an edge e in Ga not in the tree, there is a unique positively oriented 
circuit c e whose edges are either e or edges in the tree. The relation w(c e ) = 
a(c e )p^ Ce ^ q H ( Ce \ which u must satisfy from condition (d) above, now imposes the 
value of the weight of e. This defines a weight 00 on Ga- 

If e is in A, then c e is also in A and therefore 0(c e ) = 0: the weight assigned to 
e thus belongs to A*, verifying condition (a). If instead e crosses £ counterclock- 
wise (resp. clockwise), we have 0(c e ) = 1 (resp. —1), verifying in turn condition 
(b) (resp. condition (c)). 

Set W : Zi(A) -> A* to be 

W(c) =uj(c)a(c)p" i>{c) q- v{c) . 

This function is a group homomorphism and equals 1 on all positively oriented 
circuits c e , which generate Zi(A) (a fact which is left to the reader). Thus, W 
is constant equal to 1, proving that condition (d) holds for arbitrary oriented 
circuits. 

We can now infer the value of u(t) for a tiling t: given a fixed base tiling to, 
write 

u(t) = u(to)u(t - t ) = uj(t )a(t - * )/ (t -*°V ( *-* o) = uj(t )a(t; to)p^ q^ . 

Clearly, all Kasteleyn weights are obtained by this construction. ■ 

Consider a periodic track segment A with n black and n white squares, A = 
clA, £ the induced cut and a Kasteleyn weight uo on Ga- We now construct 
the associated n x n Kasteleyn matrix Ma with coefficients in A Ptq . Begin by 
separately labeling black and white squares of A (i.e., vertices of Ga)', more 
formally, let rjb (resp. i] w ) be a bijection from {1,2, ... , n} to the set of black 
(resp. white) squares. A pair r\ = (i] b , r) w ) allows us to prescribe a specific sign to 
a tiling t: a(t; rj) is the sign of the permutation 

Vw loto Vb ■ {!,••■ ,n} -> {1,... ,n}. 

In the formula above, t is interpreted as a bijection from black squares to white 
squares of A. Thus, a tiling t admits a sign once either a base tiling t or a 
labeling 77 is given: we have a(t; to) = a(t; rj)/a(to] rj). 

Set (M^)ij to be the weight of the edge joining i]b(j) and i] w (i), or zero if no 
such edge exists. From conditions (a), (b) and (c) in the definition of a Kasteleyn 
weight, we may write the Kasteleyn matrix Ma = p~ l N~^+MA+pN~^ uniquely for 
three matrices with entries in A q : (N^, Ma, N£) is the Kasteleyn triple associated 
to the Kasteleyn weight. Non-zero entries in (resp. iVj) correspond to edges 
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crossing £ clockwise (resp. counterclockwise); non-zero entries in M correspond 
to edges in A. Notice that a zig-zag cut (or more generally a wiggling cut) yields 
a null or N£. 

Let $a(p, q) = det Ma(j), q) be the q-flux polynomial of A or A. This notation 
implicitly assumes two facts: the innocuous dependence on the choice of cut and 
Kasteleyn weight. More precisely, if A — clA = clA' and Kasteleyn weights 
are assigned to A and A' then det(p~ 1 N^ f + M A > + pN£,) may be obtained 
from det(p _1 iV^ + M A + pN£) by multiplication by a factor of the form ±p*q*. 
This follows from the proposition below, which also provides a combinatorial 
interpretation for $a(p, q)- 

Proposition 4.6 Let A be a periodic track segment, A = c\A, fix a Kasteleyn 
weight 00 on Ga, a pair of bisections rj, a base tiling to G Ta and let Ma be 
the correponding Kasteleyn matrix. Let $a(p, q) = det Ma (p, q) be the q-flux 
polynomial of A. We then have, 

$ A (p,q) = a(t ; V )iu(t )J2p mt0) Q u{t '' t0) - 

ter A 

Furthermore, the special values $a(p, 1) an d &a(p, — 1) agree with the flux poly- 
nomials $a(p) and defined in the introduction up to a multiplicative factor 
of the form ±p* . 

Proof: Tilings of A correspond to non-zero monomials in the expansion of 
det Ma- Thus, for a fixed base tiling t , 



det M A = vMt) = <r(t ; vMt ) ^/^V 



(t-t ) 



the desired formula for $a(p, Substituting 1 and —1 for q yields the in- 
terpretations for $a(p, 1) and $a(p, — 1) which were called <E>^ and $^ in the 
introduction. ■ 
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Figure 8: Annulus and Kasteleyn weight 
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In the left of Figure |8] we show a periodic track segment A; the attachments 
are the vertical sides and the outer boundary is the bottom side. In the right, we 
show a Kasteleyn weight constructed as in Proposition ^4.5| by assigning weight 1 
to the edges of the maximal tree indicated by solid lines. Thick lines indicate a 
tiling t - Dotted lines indicate edges not belonging to the tree. We thus obtain a 
12 x 12 Kasteleyn matrix whose determinant is 

<Mp,g) = <r 1 V 2 (<? 3 V + 

(g 36 + 3 g 35 + 3 q u + 4g 33 + 6 g 32 + 6 g 31 + 7g 30 + 6 q 29 + 6 g 28 + 7 q 27 + 
Qq 26 + 6g 25 + 7g 24 + 6g 23 + 6g 22 + 4g 21 + 3g 20 + 3g 19 + g 18 ) p 3 + 

(g 30 + 3 g 29 + 3 g 28 + 4 g 27 + 9 g 26 + 12 g 25 + 16 g 24 + 24 g 23 + 33 g 22 + 
41 g 21 + 45 g 20 + 51 g 19 + 57 g 18 + 51 g 17 + 45 g 16 + 41 g 15 + 
33 g 14 + 24 g 13 + 16 g 12 + 12 g 11 + 9 g 10 + 4 g 9 + 3 g 8 + 3 g 7 + g 6 ) p 2 + 

(g 18 + 3g 17 + 3g 16 + 4g 15 + 6g 14 + 6g 13 + 7g 12 + 6g n + 6g 10 + 7g 9 + 
6g 8 + 6g 7 + 7g 6 + 6g 5 + 6g 4 + 4g 3 + 3g 2 + 3g + l)p+l) 

Also, $^ (p) = p~ 2 (p A + 91 p 3 + 541 p 2 + 91 p + 1) (with roots approximately equal 
to -84.619, -6.2077, -.16109, -.011818) and $ A (p) = p~ 2 (p 4 + p 3 + p 2 +p+ 1), 



in agreement with Theorems 1.1 and [L3 



5 Kasteleyn matrices for covers 

In this section we present an algebraic counterpart of the juxtaposition A™ of 
n copies of a periodic track section A, from which we obtain some preliminary 
information concerning the roots of $^(-,g). 

Proposition 5.1 Let A be a periodic track segment with a Kasteleyn weight in 
A = cl(A) with triple (N~,M,N + ) = (N±(q),M A (q),N£(q)). Then there is a 
natural Kasteleyn weight in A n = cl(A n ) with triple (N^ n , M^, N^ n ), where 
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Thus, M A n (p, g) = p- 1 N An (g) + M A ~ (?) + pN^ n (g) . 



(■■■ (-l) n+1 N~\ 
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Proof: Here the vertices in A n = Ao • • • A n _i are labeled in the following order: 
first those belonging to A (with the same order applied in the Kasteleyn triple 
for A), then those in Ai (again with the same order) and so on. The Kasteleyn 
weight in A™ assigns to each edge contained in some A& the same weight as in A. 
For edges trespassing a cut £ fc+ i clockwise (resp. counterclockwise), < k < n—1, 
we assign as weight the corresponding entry in (resp. iVj); notice that these 
weights belong to A*. Finally, for an edge in A n trespassing £_i = £ n _i clockwise 
(resp. counterclockwise), we assign as weight (— l) n+1 p _1 (resp. (— l) n+1 p) times 
the corresponding entry in (resp. N£). This clearly yields the matrices in the 
statement; we are left with proving that this construction obtains a Kasteleyn 
weight in G^n. 



Conditions (a), (b) and (c) in Definition [4.4] are clearly satisfied. We check 
(d) for the hole basis. For small holes, this is easy. Let £ be the large hole in 
Ga and its n-cover £ n be the large hole in G^n. For the Kasteleyn weight in A, 
ou(£) = a(£)p^q u ^ = (— l) k+1 p, where the boundary of the large hole has length 
2k (see Definition [p). We have a(Q = (-l) nfc+1 , 0(4) = 1 and u(£ n ) = and 
we must check that u(£ n ) = (—l) nk+1 p, which is the case by construction. ■ 

For a Laurent polynomial P(X) = a(X — Ai) • ■ • (X — \ m )X b , Ai, . . . , X m ^ 0, 
the n-th root power is 

P [n \X) = a n {X - (-l) n+1 A?) • ■ ■ {X - (-l) n+1 A™ )X b 

= n p(c 2j+l n 

j=0,l,~ >«-l 

where ( = — exp(ni/n) and Y n = X. 

The reader may check that in the product above, the terms containing Y k , 
k ^ (mod n), cancel out. Let 

^kO^l) ■ ■ ■ > J'rri) ^ 3Cii ' ' ' •Ei^ 

l<il<---<ifc<m 

be the usual symmetric functions. If P{X) = (a m X m + ■ • • + a )X b = a m (X — 
AO • • • (X - \ m )X b , then clearly pW(X) = (at ] X m + ■■■ + a [ Q ] )X b , where 

4 H = (-l) m - fc a> m - fc ((-l) n+1 A^, . . . , (-l) n+1 X n J 
= a^cr m _fc((— Ai) n , . . . , (— A m ) n ). 



Proposition 5.2 Let A be a periodic track segment, A = cl A, q a fixed complex 
number, P\{p) = $a(p, q) and P n (p) = &A n (p, <?)■ Then P n is the n-th root power 
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For most q, in particular for q = 1, the number of non-zero roots of Pi(p) is 
/max — /min (recall that / max (resp. / m i n ) is the maximum (resp. minimum) flow 
for tilings of A). However, for special values of q some of the extremal coefficients 
of Pi(p) may vanish thus causing the number of non-zero roots to go down. 

Proof: Let ( = — exp(7ci/n) and r n = p. Let M be the block diagonal matrix 
with blocks (C 2k+1 ry 1 N A (q) + M A (q) + (C 2k+1 r)NX(q), k = 0, 1, . . . ,n-l whose 
determinant is clearly p[ n \ The reader may check that M = X- 1 M A n{p,q)X 
where X is the product of the block diagonal matrix diag(J, (rl, . . . , £ n_1 r n-1 J) 
and the discrete Fourier transform block matrix, whose (j, fc)-block is ( 2 ^ h I, j, k = 
0, ... ,n- 1. ■ 

The above proposition already yields considerable information concerning the 
roots of the polynomials = &a(p, ±1)- F° r instance, if Ai, . . . , A m are the 

non-zero roots of we must then have <Tfc((— 1)™ +1 A™, . . . , (— l) n+1 A^J > for 
all k and n, since this number counts tilings with a given flux in cl(A n ) (up to a 
positive factor). In order to prove the main theorems, however, we need sharper 
estimates of such numbers; these require the tools of the next sections. 



6 Tilings of track segments and — 1-counting 

A tiling t of a track segment A is a decomposition of A as a disjoint union of 
squares with a preferred edge contained in cto or ol\ and dominoes. The preferred 
edge indicates the position of the missing half-domino. The reader will have no 
difficulty in extending the tiling of the track segment in Figure |^ to the right but 
will notice that no extension to the left exists: not every tiling of a track segment 
is a restriction of a tiling of a larger region. 



Figure 9: Tiling of a track segment, with indices 

The flux 4>(t, «o) of a tiling t of a track segment A relative to cto counts 
preferred edges in ao with sign: preferred edges which are sides of white (resp. 
black) squares count positively (resp. negatively). Similarly, <p(t, ati) counts 
preferred edges in a.\ but now black and white count positively and negatively, 
respectively. The difference <f>(t, ai) — 4>(t, cto) is the number of black squares in 
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A minus the number of white squares in A and therefore does not depend on the 



tiling t: notice the similarity with Lemma |27T . 

We now associate a height function 9 to a tiling t of a track segment A; 
as before, we choose a base vertex Vh in (3 Q and set 9{yh) = 0. Extend 9 to 
the remaining vertices of A by local instructions: it is forbidden to walk along 
preferred edges or edges which are trespassed by dominoes. The rest is similar: 9 
goes up by 1 along an allowed edge if there is a white square to its left (or a black 
square to its right); conversely, 9 goes down by 1 if there is a black square to its 
left (or a white square to its right). Two height functions for the same tiling with 
different base vertex differ by an additive integer constant. 

For a fixed base vertex, the restrictions of height functions to (3 do not depend 
on the tiling t. Choose a reference vertex v r in {3f as in Lemma |3^, 9(v r ) = 
4</>(t, ao) + c for some constant c. The value of a height function on $ therefore 
depends only on the flux of t. 



A characterization analogous to Proposition |3J] is immediate. Let 9 be an 
integer valued function defined on the set of vertices of A. Then 9 is the height 
function of some tiling t if and only if the following conditions hold: 



(a) 8(v b ) = 0; 

(b) if the oriented edge joining vertices vo to v i is in $ or [3 and there is a white 
(resp. black) square to its left then 9{vi) — 9{vq) = 1 (resp. —1); 

(c) if the oriented edge joining vertices vq to v\ is in the interior of A or in 
an attachment and there is a white (resp. black) square to its left then 
9{yi) — 9{vq) = 1 or — 3 (resp. —1 or 3). 

The volume v{t) of a tiling t of a track segment A is a weighted sum of values 
of 9(v) over the vertices of A. Vertices in j3 and $ have weight and interior 
vertices have weight 1. A vertex v in the interior of an attachment has weight n/4 
if it belongs to n squares. The volume of a track segment is thus a quarter-integer. 

Tilings t and t! of A and A' can be juxtaposed to produce a tiling of AA' if 
and only if a± and a' have the same shape and the sets of preferred edges of a± 
and a' Q coincide. The subset of preferred edges of a tiling t at a (resp. a\) is 
the 0-index (resp. 1-index) of t: thus, t and t' can be juxtaposed if the 1-index of 
t coincides with the 0-index of t! . More generally, a 0-index (resp. 1-index) is a 
subset of the set of n edges of a>\ (resp. the set of rii edges of ao). Not all indices 
arise as the index of a tiling of A: for instance, the index of a tiling cannot have 
two preferred edges which are sides of the same square. Notice that a tiling of 
AA' can be obtained in exactly one way as a juxtaposition tt' of tilings t of A 
and t' of A'. A tiling t of A can be closed-up to yield a tiling clt of clA if its 
two indices are equal. Either index of a tiling t of A determines the restriction of 



24 



Figure 10: Indices of flux 1 



the height function 9 to the corresponding attachment, and therefore the flux of 
t relative to the attachment. Finally, a tiling t of A and a tiling t' of A' can be 
juxtaposed only if <fi(t, a\) = <fi(t', a' ) (this is however not a sufficient condition). 

For the next definition, it is convenient to label indices. Order the 0-indices 
(resp. 1-indices) so that the associated fluxes are non-decreasing; among in- 
dices with the same flux, use the order induced by the natural bijection with 
{0, . . . , 2 n ° - 1} (resp. {0, . . . , 2 ni - 1}). This assigns labels from {1, . . . , 2 n °} 
and {1, . . • , 2" 1 } to the sets of indices which are non- decreasing in flux. On the 
first (resp. second) line of Figure 110 we show all 0-indices (resp. 1-indices) with 



associated flux / = 1 of a track segment (the sign of the flux corresponds to 
painting the lower right hand corner white). 

Definition 6.1 Let A be a track segment with attachments «o and a± of length 
n and n\. The 2 n ° x 2™ 1 connection matrix has its entries in C^ 1 / 4 , <7 -1 / 4 ]. 
The (io,ii) entry is the sum of q u ^ over all tilings t of A with indices labelled ig 
and %\. 



In particular, the entry (io,ii) is if the indices labelled z'o and i\ induce 
different values of </>(•, a>i) and 0(-, a' ). The matrix C*a thus splits into rectangular 
blocks Ca,/, arranged diagonally, labelled by the constant value / of the flux on 
the block. 

We provide a more concrete interpretation for the entries of the connection 
matrix. Given a pair of indices (zo,^i), we show how to prune the track segment 
A to obtain a quadriculated disk D io il . Squares of A which are not adjacent to 
attachments will be left undisturbed in D i0til . Consider now each square adjacent 
to an attachment: if no side of this square belongs to either index, the square 
also belongs to D i0iil ] if exactly one side belongs an index, the square is not 
included in D ioil ; finally, if two or more sides belong to indices, D i0jil is not 
defined. Clearly, when D io>il is defined, tilings of D io>il are in natural bijection 
with tilings of A with indices i$ and i±; also, when D i(itil is undefined, there are 
no tilings of A with indices io and i\. The converse is not true: is defined 
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even when the corresponding fluxes are different but will not be balanced and 
will therefore admit no tilings. Summing up, the (io,ii) entry of Ca is, up to 
multiplication by an integer power of g 1//4 , the g-counting of tilings of D i0jil and 
if Di (h i 1 is undefined. 

Strictly speaking, the matrix Ca depends on the choice of the base vertex. 
Different choices, however, merely have the effect of multiplying the connection 
matrix by a power of g 1 / 4 . 

Connection matrices behave well with respect to closing-up and juxtaposition. 

Proposition 6.2 Let A be a periodic track segment (with prescribed base and 
reference vertices). Then §a{p, q) — P*Q*YlfP^ r C'Aj, where $a is the q-flux 
polynomial of the annulus cl A and p* and q* denote arbitrary integer powers of 
p 1 / 4 and g 1 / 4 . 

Let A' and A" be track segments (with prescribed base and reference ver- 
tices) such that the attachments a[ and have the same shape. Then Ca>a" = 
q*CA'CA", where the track segment A' A" receives arbitrary base and reference 
vertices and q* denotes an arbitrary integer power of g 1//4 . 

Proof: For an appropriate choice of base and reference points V} } and v r , a 
tiling c\t of clA induces a tiling t of A with u(t;v b ) = u(c\t;Vb), 4>(t; i>&, v r ) = 
(f)(c\t;Vb,v r ) and 0- and 1-indices of t are equal. Conversely, a tiling t of A with 
equal 0- and 1-indices induces a tiling clt of cl A. 

For the second statement, we may assume the base and reference vertices 
to belong to the common attachment a[ = ckq since this can only modify the 
final result by multiplying by an integer power of g 1//4 . Notice that with this 
choice of base and reference vertices, the set of values of the flux / is the same 
for both track segments. The volume of a tiling t't" of A' A" clearly satisfies 
v{t't") = u(t') + v(t"). The result now follows from the definition of matrix 
product. ■ 

Corollary 6.3 Let A be a periodic track segment (with prescribed base and ref- 
erence vertices), q a fixed complex number and Ai, . . . , A m the non-zero roots of 

& A (p, q) = ap /min (p m + a m -iP m_1 • • • + ao) 
= ap f ^{p-X 1 )---{p-X m ). 

Then 

trCS i/ = aV /M _ / ((-A 1 r,...,(-A m r). 

We are ready for the main two results of the section. The matrix j is 
obtained from Ca,/ by substituting q — — 1. 
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Theorem 6.4 Let A be a periodic track segment. The eigenvalues of the square 
matrices * are or roots of unity. Furthermore, for any non-zero eigenvalue 
the algebraic and geometric multiplicity are equal. 

The algebraic multiplicity of the eigenvalue A = (i.e., the multiplicity of 
the root zero in the characteristic polynomial of j) may be different from its 
geometric multiplicity (i.e., from dimker((7^ j — XI))- For the example in Figure 
PH| (for which / = 1) we take as base point the central point on the top (/3 ) and 
we then have 

/0 0\ 
10 1 
10 1 
\0 0/ 

which is nilpotent. Nonzero entries correspond to classes of tilings with a single 
element of even integer height. Some zero entries correspond to empty classes 
and other correspond to classes where cancellation occurs. 

Proof: Let C = ^ so that C n = C^ n j. The entry io,ii of C n , n 6 
{1,2,3,...}, equals, up to multiplication by an integer power of (— 1) 1//4 , or 
the — 1-counting of tilings of the quadriculated disk (A n ) i0jil . By Theorem |1.2j in 
the introduction, these coefficients can only assume a finite number of different 
values. Thus, for some distinct natural numbers a and b we have C a = C b , from 
which the result follows. ■ 
We are now ready to prove Theorem |1.3|. 



Theorem 6.5 Let A be a balanced quadriculated annulus. All non-zero roots of 
the signed flux polynomial $^(p) are roots of unity. 

Proof: Call these roots Ai,...A m . Let A be a periodic track segment with 
A = cl A. From Theorem ^]4| we have that, for any fixed k, the sequence 

^((-Ai) B , • • • , (-A m ) n ) 

is eventually periodic in the variable n. In other words, there exist positive 
integers n Q and n\ such that 

^((-Ai)^ ftl , . . . , (-X m ) n+Nni ) = o-kd-X!)* 1 , . . . , (-X m ) n ) 

for all k, all non-negative iV and all n > n . In particular, 

^(((-AO" " 1 ) 2 , . . . , ((-A m )"^) 2 ) = ^(((-AO^), . . . , ((-AJ" " 1 )) 

for all k, meaning that taking squares effects a permutation of multiset 

{(-Xi) noni , • • • , (-X m ) noni }. 

The ml power of this permutation is the identity, implying that, for all j, 

(-Xj) 2m '" n ° ni = (-A J ) nori1 

and thus Xj is a root of unity. ■ 
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The reader may wonder what would have happened if instead of A we had 
chosen A' with cl A = cl A' = A. Actually, arbitrary powers of Ca and Ca' have 
the same trace since both traces count the tilings of A n . More directly, as the 
reader can easily convince himself with a picture, for sufficiently large n there are 
(non-periodic) track segments Ai and A2 with A1A2 = A n , A2A1 = (A')". 

The estimate on the algebraic degree of Xj produced in the above proof is 
ridiculously large: Xj, being a root of an integer polynomial of degree at most 
/max - /min has algebraic degree at most / max - / mm . 



7 Eventual compatibility and ^-counting, q > 

Two tilings t\ and £2 of a quadriculated surface are said to differ by a flip if 
they are identical except at a 2 x 2 square, where one tiling has two horizontal 
dominoes and the other has two vertical dominoes. Two tilings ti and t 2 of a 
disk or band differ by a flip if and only if their height functions 9\ and 9 2 differ at 
a single point, the center v of the square, and then 9 2 {y) — Oi(v) = ±4. Assume 
9 2 {v) — 9\{v) = 4: from the local characterization of height functions, v is a 
local minimum (resp. maximum) of 9i (resp. 9 2 ). Moreover, flips are admissible 
exactly at local extrema of a height function. 

It follows from these remarks ( |TD| , |J ) that any two tilings in a quadriculated 
disk can be joined by a finite sequence of flips. The situation in quadriculated 
annuli is not so simple: we state a special case of Theorem 4.1 in ||. 

Theorem 7.1 Let A be a wall-free quadriculated annulus and let t\ and t 2 be 
tilings with the same flux. Then there exists a finite sequence of flips joining t\ 
and t 2 . 



The hypothesis of A being wall-free is necessary as mentioned in ||, page 224; 
the counterexample is also given in Figure |] above. 

Let A = clA be a quadriculated annulus with universal cover A°° = 
. . . A_ 1 A A 1 .... Let ti and t 2 be tilings of A inducing periodic tilings n^fti) 
and 7r~ l (t 2 ) of A°°. A tiling t 12 of A°° interpolates from ti to t 2 if there exists an 
integer K such that t± t 2 an d 7r -1 (£i) (resp. 7r _1 (t 2 )) coincide in A^ for k < —K 
(resp. k > K). Two tilings t\ and t 2 of A are eventually compatible if they admit 
interpolations £12 and £21- 



From Lemma 2.1, tilings with different fluxes are not eventually compatible. 



Also, the tilings ti and t 2 in Figure |] have the same flux. Figure |TT] shows an 
interpolation £ 12 but no interpolation t 2 ,i exists (as is clear from the picture) and 
hence t\ and t 2 are not eventually compatible. The following proposition shows 
that this phenomenon does not occur in wall-free annuli. 
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Figure 11: Possible and impossible interpolations 



Proposition 7.2 Let A be a balanced wall-free quadriculated annulus and let t\ 
and t2 be tilings with the same flux relative to the same cut. Then t\ and ti are 
eventually compatible. 

Proof: Compatibility is an equivalence relation. Also, tilings of A differing by a 
flip are eventually compatible (see Figure [12]). The proposition now follows from 
Theorem |7.1| . ■ 



Figure 12: Tilings differing by a flip are eventually compatible 

Let A be a periodic track segment and let i be an index: i is left-active (resp. 
right-active) if there exists a tiling of the left (resp. right) half-band . . . A_iAq 
(resp. AiA 2 . . . ) with index i at £i. 

An alternative definition is the following. For a value / of the flux, we con- 
struct the f -index graph If, an oriented graph whose vertices are the indices of 
flux / at £i and we join i and i' by an oriented edge if there exists a tiling of A 
with index % at £_i and i' at £i. An index % is left-active (resp. right-active) if 
there exists an infinite path in If ending (resp. starting) at i. 

Proposition 7.3 Let A be a periodic track segment, f a value of the flux and If 
the f -index graph. Then there exists an integer N such that for all n > N and 
all indices i and i' the following condition holds: there is a path in If of length n 
starting at i and ending at i' if and only if i is right-active and %' is left-active. 

Proof: Let I be the number of vertices in If (i.e., the number of indices of flux 
/). If there is a path of length 7 + 1 starting at % then the path must pass twice 
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through some vertex and it is straightforward to construct an eventually periodic 
infinite path starting from % (an infinite path i Q ii ... is eventually periodic if there 
exist positive integers k and J such that ij + k = ij for all j > J; k is the period) . 
Thus, if a path starts at i, ends at %' and has length greater that I then i is 
right-active and i! is left-active. 

A tiling of cl(A n ) induces a circuit of length n in Z/; conversely, any such 
circuit is induced by some tiling. Let i be a right-active and i' be a left-active 
index. Construct an eventually periodic path starting at i. Let k be its period: 
the path eventually correspond to a tiling of cl(A fc ) which is compatible with 
tf,mm by Proposition [T2| (here i/, m in £ is the periodic tiling corresponding to 
the minimal height function 0/ jm i n with flux / on A°°). This shows us how to 
modify the path to obtain an eventually periodic path of period 1, corresponding 
to a tiling on the right half-band eventually coinciding with £/ jm i n . The same 
construction applied to the path ending at H yields a tiling of the left half-band 
also eventually coinciding with t/ jm i n . It is now clear how to juxtapose these two 
tilings to obtain a tiling of A . . . A n with indices i and i' for any sufficiently large 
n. ■ 

An index i is bi-active if both left- and right-active. The bi-active submatrix 
C*,a,/ is obtained by removing from Caj all rows and colums associated with 
indices which are not bi-active. 

Proposition 7.4 For any fixed complex number q the spectra o/C^a,/ and Caj 
differ by null eigenvalues. 

Proof: An index i is strictly left-active if left- but not right-active, strictly right- 
active if right- but not left-active, inactive if neither left- nor right-active. We 
permute lines and columns of C — Caj so that indices appear in the following 
order: left-active, bi-active, right-active, inactive, thus turning C in a 4 x 4 block 
matrix (with blocks of different sizes). The block 6*2,2 is C*,A,/- 

We claim that (C n ) 2 ,2 = (^2,2)" (i-e., the 2, 2 block in the n-th power of C is 
the n-th power of its 2, 2 block). Indeed, a monomial in the expansion of an entry 
of (C n )2,2 is a path in Z/ starting at a bi-active index and ending at a bi-active 
index: by definition of bi-activity this path can be extended both to the left and 
right and therefore, again by definition, all indices in the path are bi-active. In 
other words, the only non-zero monomials in the expansion of (C n )2.2 are those 



in (C 2i2 ) n . 



From Proposition |7.3|, for all n > N, 
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It is now clear that the spectra of C n and (C n )2,2 differ by null eigenvalues. From 
the previous paragraph, the spectra of C n and (C 2 ,2)™ differ by null eigenvalues 
and we are done. ■ 

A matrix X is positive if all entries are real positive. A square matrix X is 
eventually positive if X n is positive for some positive n. 

We use the following simple modification of the classical Perron theorem (f|). 

Theorem 7.5 If M is an eventually positive matrix then it admits a simple pos- 
itive eigenvalue X strictly larger than the absolute value of any other eigenvalue. 



Theorem 7.6 Let A be a periodic track segment with a wall-free closing-up 
A = cl A. For any fixed positive real number q, the square matrix C&j has a 
simple positive eigenvalue X strictly larger than the absolute value of any other 
eigenvalue. 

Proof: From Proposition |7.3| , the bi-active submatrix C* a,/ is eventually pos- 
itive. From Theorem |7.5| , C^a,/ admits a simple positive eigenvalue A strictly 
larger than the absolute value of any other eigenvalue. The result now follows 
from Proposition V 



We are ready to restate and prove Theorem 1.4. 



Theorem 7.7 Let A be a balanced wall-free quadriculated annulus and q a fixed 
positive real number. The non-zero roots (in p) of the flux polynomial ^a{p,q) 
are distinct real negative numbers. 

Proof: Call Ai, . . . A m the non-zero roots of $a(p, <?), with |Ai| > |A 2 | > ... > 
|A m | > 0. Assume by induction that Ai, . . . , A^-i are real negative and that 

|Ai| > ■ • • > |Afc_i| > |Afc| = ■ • • = |Afc'| > |Afc/ + i|, 

with k' > k. We must prove that k = k' and that A& is real negative. 



From Corollary |6.3| and Theorem 7.6 we have that the sequence of k- 



symmetric functions of powers of Xj satisfies cr fc ((— Ai) n , . . . , (— X m ) n ) = (1 + 
o(l))B%. In the notation of Corollary |6~3| , the positive constant is the 
largest eigenvalue of CA,/ max -fe (for our fixed value of q) divided by a. The ex- 
pression o~k{{— Ai) n , . . . , (— A m ) n ) is the sum of k! — k + 1 terms of the form 
((— l) fc AiA 2 • • ■ Xk-iXi) n , £ — k, . . . , k' and other terms which grow at exponen- 
tially smaller rates. Thus 

(-1)^^ • • • XU(K + ■■■ + A£) = (1 + o{l))Bl 

whence 

bl + ... + bl = l + o(l) 



31 



where be = —\e/\\e\. Thus bk, ■ ■ ■ ,bk> all belong to the unit circle and we can 
take arbitrarily large n such that 2(k' — k + 1)|&™ — 1| < 1 for all i and therefore 

|( fc '_ fc + l)_(63f + ... + 65»)|<l/ 2| 

a contradiction unless k = k'. Finally, bk — 1 and Xk is real negative. ■ 

We finish this section by estabilishing a few simple consequences of the pre- 
vious theorem. 

A sequence of real numbers is convex (resp. concave) if ak-i — 2a^ + a^+i 
is non-negative (resp. non-positive) for all k. A sequence of non-negative 
real numbers is log-concave if a\ > a k _ia k+ i for all k. In particular, log-concave 
sequences are either monotone or unimodal. 

Proposition 7.8 If the roots of P(x) = a n x n + ■ ■ ■ + a are all real negative and 
a n > then the sequence a k is log-concave. 



Proof: We proceed by induction on the degree n (the case n = 1 is trivial). The 
result now follows from checking that if the coefficients of Q(x) form a log-concave 
sequence then the coefficients of (x + c)Q(x), c > are likewise log-concave. ■ 

The converse is not true: the coefficients of X 2 + X + 1 form a log-concave 
sequence but the zeroes are not real. 

Corollary 7.9 Let q be a fixed positive real number and let aj be the coefficient 
of pf in <&ji(p,q). Then the sequence a/ is log-concave. 



Proof: This follows directly from Theorem 1.4 and the previous proposition. 



Corollary 7.10 Let bf (resp. Cf) be the highest (resp. lowest) exponent of q in 
non-zero terms of $>a(p,q) of the form pf q* . Then the sequence bf (resp. Cf) is 
concave (resp. convex). 



Proof: This follows from the previous corollary in the limit cases q — > and 
q — > +oo. ■ 
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